1 



The role of Stewartson and Ekman layers in 
turbulent rotating Rayleigh-Benard 
convection 

Rudie P.J. Kunnen 1 , Richard J. A.M. Stevens 2 , Jim Overkamp 1 , Chao 
CN Sun 2 , GertJan F. van Heijst 1 , Herman J.H. Clercx 1 ' 3 

Oh 1 Department of Physics and J.M. Burgers Centre for Fluid Dynamics, Eindhoven University of 

CD Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands 

^0 2 Department of Science and Technology and J.M. Burgers Centre for Fluid Dynamics, 

University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands, 



O 



department of Applied Mathematics, University of Twente, P.O. Box 217, 7500 AE 

Enschede, The Netherlands 



When the classical Rayleigh-Benard (RB) system is rotated about its vertical axis roughly 
three regimes can be identified. In regime I (weak rotation) the large scale circulation 
(LSC) is the dominant feature of the flow. In regime II (moderate rotation) the LSC 
is replaced by vertically aligned vortices. Regime III (strong rotation) is characterized 
by suppression of the vertical velocity fluctuations. Using results from experiments and 
direct numerical simulations of RB convection for a cell with a diameter-to-height aspect 

•^h ratio equal to one at Ra ~ 10 8 — 10 9 (Pr = 4 — 6) and < l/Ro < 25 we identified 

the characteristics of the azimuthal temperature profiles at the sidewall in the different 
regimes. In regime I the azimuthal wall temperature profile shows a cosine shape and a 
Oh vertical temperature gradient due to plumes that travel with the LSC close to the sidewall. 

In regime II and III this cosine profile disappears, but the vertical wall temperature 

— . gradient is still observed. It turns out that the vertical wall temperature gradient in 

regimes II and III has a different origin than that observed in regime I. It is caused by 
boundary layer dynamics characteristic for rotating flows, which drives a secondary flow 
that transports hot fluid up the sidewall in the lower part of the container and cold fluid 

\^ downwards along the sidewall in the top part. 

ON 

O 

1. Introduction 

• i— i The classical system to study buoyancy driven flows is the Rayleigh-Benard (RB) sys- 

tem (Ahlers et al. 2009 Lohse & Xia 2010). In this system a layer of fluid is confined 



between two horizontal plates and is heated from below and cooled from above. The RB 
system is also very suitable to study the influence of rotation on heat transport mech- 
anisms. In order to do this the RB system is rotated at an angular speed fl about its 
vertical axis. Studies about the influence of rotation on heat transport are relevant to 
understand many geophysical and astrophysical flow phenomena such as the global thcr- 
mohaline circulation, convection in the interior of gaseous planets, and convection in the 
outer layer of the Sun. Furthermore, the studies are also very relevant for optimization of 
industrial applications. Therefore, turbulent rotating convection has been studied exten- 
sively in laboratory experiments (for e xample, |Rossby||1969| |Boubnov fe Golitsyn||1986| 



Fernando et al. 1991; Zh ong et al. [19931 |L~iu & Eckc 1997; Sakai 1997 Hart & Qlsc n|1999 



Vorobieff fe E ckc 2002} |Kunnen et al . 2008 b\ Liu fc Eck c 2009; Nie mela et alpOl O; Zhong 
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Figur e 1. The scaled heat transfer Nu(l/R6)/Nu(0), with iVu the Nusselt number defined 
in Eq. ( |f .2\ , as function of 1/Ro on a logarithmic scale. Experimental and numerical data for 
Ra = 2.73 x 10 8 and Pr = 6.26 are i ndicated by red dots and open squares, respecti vely. Data 



taken f rom | Zhong et al. ( 2009[ ) and |Stevens et al.\ ( |2009[ ) and are called run E2 in |Zhong fc| 



lAhlerslpOlOp . The transition between the different regimes, see text, is indicated by the vertical 



dashed lines. 



fc Ahlers|2010 ) , direct numerical simulations (for exam ple, |Julien et aL|19966|a Kunnen 



et aLI2006[|Schmitz fe Tilgner|2009l |Kunnen et aZ.|2009||Schmitz fe Tilgner|20io[|Stevens 
et aL|| 2010c|«[), and with combined experimental and numerical investigati ons (|Kunnen 



et al. 



2010). 



2008aHZhong et aZ.||2009l |Stevens aT1[2009| |Kunnen et aT1|2010a|6[ |Weiss et al. 



In this paper we indicate the rotation rate around the vertical axis by the Rossby 
number Ro, which compares the inertial and Coriolis forces in the system. Here, Ro is 
defined as 



Ro 



Ra 



1 



(1.1) 



PrTa 

where Pr — v/k is the Prandtl number, Ta = (2flL 2 /iy) 2 is the Taylor number and 
Ra = f3gL 3 A/(vK) is the Rayleigh number, with the rotation rate, j3 the thermal 
expansion coefficient, g the gravitational acceleration, v the kinematic viscosity, k the 
thermal diffusion coefficient, and A the temperature difference between the two plates. 
Throughout the paper we will consider a cell with aspect ratio T = D/L = 1, with L 
the height of the RB cell and D its diameter. The dimensionless heat transport in the 
system is indicated by the Nusselt number, 



Nu 



QL 
AA 



(1.2) 



where Q is the heat-current density and A the thermal conductivity of the fluid in the 
absence of convection. 

In figure [T] we show a typical measurement of the heat transport enhancement with re- 
spect to the non-rotating case as function of the rotation rate ( Zhong et aZ.|2009 Stevens 
et al. 2009). The figure shows that, depending on the rotation rate, three different regimes 
can be identified. Regime I (weak rotation), where no heat transport enhancement is ob- 
served, regime II (moderate rotation), where a strong heat transport enhancement is 
found, and regime III (strong rotation), where the heat transport starts to decrease. We 
note that the division between regime I and regime II is obvious as there is a bifurca- 



tion (Zhong et al. 2009| Stevens et al. 2009 Weiss et al. 2010). Furthermore, we note 



that the log plot in figure [1] makes the transition from regime II to regime III look more 
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sudden than it actually is. Zhong et al. (2009) and Stevens et al. (2009) have shown that 



the position of the onset does not depend on Ra and Pr. However, the position of the 
maximum heat transfer enhancement shifts towards lower 1/Ro when Ra is increased 
or Pr is decreased. In addition, the maximum heat transfer enhancement decreases with 
increasing Ra and decreasing Pr. For a detailed discussion on the influence of Ra and Pr 



we refer the reader to Zhong et al. (2009) and Stevens et al. (2010c). Flow visualization 



experiments (Boubnov & Golitsyn 1990 Kunnen et al. 2010a) and the analysis of the 



flow structures obtained in numerical simulations (Stevens et al. 2009 Kunnen et al. 



2010 a) have confirmed that this division of regimes coincides with changes observed in 
the flow patterns and flow characteristics. 

In regime I (1/Ro < 0.5) the large scale circulation (LSC), typical for non-rotating 
RB convection, is still present, because the Coriolis force is too weak to overcome the 
buoyancy force that causes the LSC. Stevens et al. ( 2009[ ) showed that there is a sharp 
transition to the regime where rotational effects become important, while at the same 



time the strength of the LSC is decreasing (Kunnen et al. 2008a). Zhong & Ahlers 



(2010) experimentally found that the time-averaged LSC amplitudes decrease strongly 
at the transition from regime I to regime II, see figure 13 of their paper. In regime II 
(0.5 < 1/Ro < 6.67) the LSC is replaced by predominantly vertically oriented vortical 
columns as the dominant flow structures and a large increase in the heat transport is 
observed. 



This enhanced heat transport has been ascribed to Ekman pumping 


Rossby 1969 


Julien et aL|1996a||Vorobieff & Ecke|2002||Kunnen et aZ.|2008a||King et al. 


2009 Stevens 


et aL|2009 1. The effect of Ekman pumping, and thus the observed heat transfer enhance- 


ment, depends strongly on the Ra- and Pr- values (Zhong et al. 2009 


Stevens et al. 



2010c). When the rotation rate is increased further a large decrease in the heat trans- 
port is observed, because the vertical velocity fluctuations are suppressed due to the 
rotation ( Kunnen et aL|[2~0086 2010a I. We will call this regime III. 

In this paper we address the question of how these different regimes can be identi- 
fied from measurements of the azimuthal wall temperature distribution and the vertical 
temperature gradient along the sidewall, with probes that are embedded in the sidewall 
of a RB convection cell. This method of sidewall temperature measurements has been 



introduced by Brown et al. ( 20056[) and a validatio n of this method is described in the 



second paragraph of Section 2 of Brown & Ahlers (20066). Recently, Zhong & Ahlers 



(2010) have extensively studied the properties of the LSC in RRB convection in aspect 
ratio r = 1 experiments using this method. These measurements covered the Ra number 
range 3 x 10 s < Ra < 2 x 10 10 , the Pr number range 3.0 < Pr < 6.4, and the 1/Ro 
number range < 1/Ro < 20. In these measurements detailed statistics about the ther- 
mal LSC amplitude (i.e. the amplitude of the cosine fit to the azimuthal temperature 
profile at the sidewall), the LSC orientation over time, the temperature gradient along 
the sidewall, the retrograde rotation of the LSC, the frequency of cessations, etc., over 
this wide parameter range were determined. 

This paper is organized as follows. The experimental setup that has been built at the 
Fluid Dynamics Laboratory at Eindhoven University of Technology is discussed in sec- 



tion 2. The design of this setup is closely based on the Santa Barbara design (Brown 
et al. 2005a). Subsequently we will present some experimental results on the properties 
of the azimuthal temperature profiles in section 3. For a much more detailed set of ex- 
perimental results on the properties of the azimuthal temperature profiles we refer to the 



paper of Zhong & Ahlers (2010). In Section 4 we compare the experimental data with the 
azimuthal temperature and vertical-velocity profiles close to the sidewall found in direct 
numerical simulations (DNS). In order to explain the sidewall temperature measurements 
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Figure 2. Schematic diagram (not to scale) of the Eindhoven RB apparatus. From bottom 
to top the figure shows the support plate (A), an insulation layer (B), the bottom adiabatic 
shield (C), an insulation layer (D), the bottom copper plate (E), the Plexiglas sidewall (F), the 
adiabatic side shield (G), the top copper plate (H), and the Plexiglas top plate (I). See further 
details in the text. 



we study the azimuthally averaged flow profiles obtained from DNS. The analysis pre- 
sented in section 5 reveals that the sidewall temperature measurements, particularly the 



presence of the vertical wall temperature gradients in regimes II and III ( Zhong & Ahlers 



2010), can be explained by the presence of Stewartson layers that are formed along the 
sidewall. 



2. Experimental setup 

Based on the setup described by Brown et al. (2005a) we have built a new RB setup 



in Eindhoven, which is suitable for high precision heat transport measurements. The 
convection cell has a diameter D and height L of 250 mm, making the aspect ratio 
T = D/L = 1.0. The modular design of the setup provides the ability to perform mea- 
surements at different aspect ratios as well, although this paper will only describe the 
r = 1 setup in detail. During rotating experiments the RB cell is placed on the Eindhoven 
Rotating Table Facility (RTF) (for details, see van Bokhovcn 2007). A rotating connec- 
tion is available for coolant fluid, with separated in and out flow tubes. The rotation rate 
fl of the RTF can be controlled from to 10 rad/s with a resolution of 0.001 rad/s and 
a relative accuracy of 0.5%. The rotation rates were kept below f2 = 1.57 rad/s, yielding 
Froude numbers Fr — il 2 (D/2)/g less than 0.03 for all runs and much smaller for most. 
This means that the effect of centrifugal forces can be considered negligible. 

A schematic diagram of the RB cell is shown in figure [2] From bottom to top, we first 
find the support plate A (400 x 400 mm aluminium), used to mount the setup on the 
RTF and to align the rotational axes of the table and the convection cell so that the 
offset of the cylinder axis and the rotation axis is less than 0.1 mm. Next, part B is a 
10 mm insulation layer to prevent heat loss of the bottom adiabatic shield C (310 mm 
outer diameter, 10 mm thick aluminium) to the support plate. The shield is fitted with a 
250 W heater to actively control its temperature. Inside this adiabatic shield, part D is 
another 10 mm insulation layer, followed by the copper bottom plate E (270 mm outer 
diameter, 30 mm thick) of the convection cell. The back of the bottom plate is covered 
uniformly with two double spiral grooves of 2 mm depth, 2 mm width and 6 mm spacing. 
Two 4 m, 12.3 fl resistance wires are epoxied into these grooves. From the back of the 
plate, five small holes (one at the centre, two at a radius of 98 mm and two at 100 mm, 
i.e. in the space that is available between the heater wires that are placed inside the 
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plate) are drilled to within 0.7 mm of its top surface and thermistors are mounted in 
these holes. 

On top of the bottom plate, the Plexiglas sidewall F with an inner diameter of 250 mm 
is placed. The thickness of the Plexiglas sidewall is 10 mm to provide enough strength 
in strongly RRB experiments. A rubber O-ring seals the interface between the bottom 
plate and the sidewall and prevents any leakage of fluid. The sidewall contains a fluid 
inlet near the bottom plate and, at opposite angular position, an outlet adjacent to 
the top plate. Care has been taken to avoid the entrapment of air in the convection 
cell, both when filling the cell and when heating the fluid. The sidewall is surrounded 
by an adiabatic side shield G made of two 3 mm thick aluminium plates. The shield 
is actively temperature controlled, ensuring that its temperature is always close to the 
mean temperature of the system. Additional insulation is present between the Plexiglas 
side wall and the aluminium side shield. 

The top of the convection cell is formed by the copper top plate H, which has similar 
dimensions as the bottom plate. The top plate contains a double spiral water channel of 
8 mm width, 26 mm depth and 25 mm spacing. The water, coming from a refrigerated 
circulator (Thermo Scientific HAAKE DC50-K41) at 12.5 L/min, cools the top plate 
down to the desired temperature. From the top of the plate, five small holes (one at 
the centre, four at a radius of 100 mm) were drilled to within 0.7 mm of the copper- 
fluid interface and thermistors are mounted in these holes. The top plate is covered by 
a Plexiglas top / and the water channel is sealed with a rubber ring. The visual access 
provides a way to check the proper operation of the cooling system. An aluminium 
construction ring is connected to the upper side of the top plate. It is supported by 6 
stainless steel support poles. With these poles, the RB cell can be fixed to achieve a fully 
watertight connection between the plates and the sidewall. To decrease the effects of air 
convection near the convection cell, the entire construction is placed in a wooden box 
covered with a 40 mm insulating layer on the inside. 

The apparatus contains 40 thermistors, which were calibrated simultaneously in a 
separate apparatus against a laboratory standard. During the calibration procedure all 
thermistors and the laboratory standard are placed in an extremely well temperature- 
controlled batch (temperature differences less than 0.002 K) to determine the resistance 
at a set of known temperatures. From this calibration data fits were made with 5-th 
order logarithmic polynomials to calculate temperature values from 3-wire resistance 
measurements. Deviations from these fits are generally less than 0.001 K. During the 
calibration procedure all thermistors are connected to the same software and hardware as 
in the real experiments. 24 of these thermistors are placed in the sidewall, forming 3 rings 
of 8 equally spaced sensors at heights 0.25L, 0.5L and 0.75L. Both the top and bottom 
plate contain 5 sensors to monitor the temperature of the plate. The remaining sensors 
are used to control the temperature of the bottom adiabatic shield and of the insulation 
around the Plexiglas sidewall. Readings of all thermometer resistances and of the bottom 
plate heater current and heater voltage were taken every second. The top (bottom) 
temperature T t (Tf,) was set equal to the area average of the five thermometers embedded 
in the top (bottom) plate. For any given data point, measurements over typically the first 
four hours were discarded to avoid transients, and data taken over an additional period 
of at least another eight hours were averaged to get the heat-current density Q, and the 
temperatures and T t . Because of the imperfect temperature uniformity of the bottom 
shield there was a small parasitic heat loss from the bottom plate of about 0.1 W, which 
was determined by measuring the required power to keep the bottom and top plate at 
40°C, a measurement which takes two days. For each measurement this parasitic heat 
loss was subtracted from the measured power. 
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Figure 3. The ratio of the Nusselt number Nu(l/Ro) in the presence of rotation to Nu(0) for 
Pr = 4.38 (T m = 40.00°C) and F = 1. The solid circles are the data obtained in the Ein dhoven 
RB setup and the open squares are the corresponding data obtained in Santa Barbara (fZhongl 



et aZ.]|2009 



Stevens et al.||2009| ). Black solid circles: Ra = 2.99 x 10 8 (A = 0.50 K). Red solid 



open squares: Ra = 5.6 x 10 8 (A = 1.00 K, run E4 in 


Zhong & Ahlers|2010 


Ra = 1.2 x 10 9 (A = 2.00 K, run E5 in 


Zhong & Ah 


ers|2010i. 



(A = 2.00 K). Red 



The measurements have been verified against the experiments of Funfschilling et ai\ 
(2005), Ahlers et al. (2009) and Stevens et ( [2009 ). For the present measurements the 
difference is in the order of 1%, well within the differences observed between different 
RB setups ( |Ahlers et aT\ \2009). In figure [3] we show the heat transfer enhancement 
Nu(Q)/Nu(0) as function of the rotation rate 1/Ro for Ra = 2.99 x 10 8 , Ra = 5.88 x 10 8 
and Ra — 1.16 x 10 9 at Pr — 4.38, which corresponds with a mean temperature of the 
fluid of 40.00°C, and compare them with similar heat transport measurements reported 
in the literature ( Zhong et aZ.||2009 Stevens et aL||2009 ). We find that the two datasets 
involving rotation also agree within 1%, which can be considered a very good agreement. 



3. Experimental sidewall measurements 

In RB experiments it is common to determine the flow properties by analyzing the 
azimuthal wall temperature profile obtained by thermistors that are embedded in the 



sidewall (Ahlers et al. 2009) and in our setup we do the same. Following Stevens et al. 



(20106) we define the relative LSC strength at midheight (S m ), based on the energy in 
the different modes of the azimuthal temperature profile, as 




/ 1- 



.0 



(3.1) 



Here E\ indicates the sum of the energy in the first Fourier mode over time, i.e. 

from the beginning of the simulation t = % to the end of the simulation t — t e , y^ t ° E to t 
the sum of the total energy in all Fourier modes over time, and N the total number 
of Fourier modes that can be determined. The relative LSC strength S m always has a 
value between and 1; here 1 indicates that the azimuthal profile is a pure cosine profile, 
which is a signature of the LSC according to Brown & Ahlers (2006&), and indicates 
that the magnitude of the cosine mode is equal to (or weaker than) the value expected 
from a random noise signal. Hence S m 0.5 indicates that a cosine fit on average is a 
reasonable approximation of the data, as then most energy in the signal resides in the 
first Fourier mode. In contrast, S m <§C 0.5 indicates that most energy resides in the higher 
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Fourier modes. Hence, we consider the LSC as dominant once S m 3> 0.5 at midheight. A 
small value of S m indicates that no single LSC is found, implying the existence of either 
multiple rolls or, in the case of rotating Rayleigh-Benard (RRB) convection, vertically 
aligned vortices. 

Figure |4](a) shows the magnitude of the different Fourier modes of the azimuthal wall 
temperature profile for the experiments at Ra = 1.16 x 10 9 with Pr = 4.38 based on the 
data of eight equally spaced thermistors placed inside the sidewall (located at z = 0.5L). 
The corresponding relative LSC strength is given in panel b. The figure clearly shows that 
the relative LSC strength is large in regime I, which indicates the presence of a LSC. 
However, for higher rotation rates, i.e. regime II, the relative LSC strength decreases 
because vertically aligned vortices become the dominant feature of the flow. For this flow 
one expects a random azimuthal temperature profile at midheight, which is confirmed 
by the low relative LSC strength S m in regime II and III, see figure Qb) . 

Extensive sidewall temperature gradient measurements for non-rotating RB convection 
were done by Brown & Ahlers (20071 and for RRB convection by Zhong & Ahlers (2010), 
see their figures 10 and 11. When the LSC is the dominant feature of the flow the vertical 
temperature gradient at the sidewall is mainly due to plumes that travel close to the 
sidewall with the LSC. Thus one would expect that the temperature gradient along the 
sidewall should decrease in regime II as the LSC disappears there. However figure 10 
of Zhong & Ahlers (20101 shows that the temperature gradient along the sidewall even 
increases in regime II and III. This nonzero temperature gradient along the sidewall in 
regime II and III is caused by the secondary flow that will be discussed in section [5j 
Further interesting information about the LSC can be found in figure 13 of |Zhong fc| 
Ahlers (2010). The figure shows that the temperature amplitude of the LSC, i.e. the 
strength of the first Fourier mode, starts to decrease at 1/Ro m 0.41 and that the LSC 
ceases to contribute significantly for 1/Ro > 0.8. Thus our experimental results in figure|4] 
are in good agreement with the results of Zhong & Ahlers ( |2010 ). It is interesting to note 
that the measurements of Zhong & Ahlers (2010) reveal an increase in the LSC strength 
between 1/Ro = and 1/Ro « 0.41. In our data this effect is not visible, since we do not 
have enough measurement points in regime I. As is discussed by Zhong & Ahlers (2010) 
there is no well-accepted theoretical explanation for this phenomenon at the moment. 



Numerical study of temperature and vertical- velocity profiles close 
to the sidewall 



We shall now consider the results of DNS of RRB convection performed at various Ro 
-- 2.73 x 10 8 and Pr = 6.26 ( see details 
10 9 and Pr = 6.4 (see details in Kunnen et al. 



values for Ra - 
for Ra = 1.00 x 



Zhong et al. 2009) and 



2010a 



which are in the 



same parameter regime as the experiments discussed above. In the simulations we solved 
the three-dimensional Navier-Stokes equations within the Boussinesq approximation in a 
cylindrical cell with T = 1 with no-slip boundary conditions at all walls, a uniform tem- 
perature at the horizontal plates, and an adiabatic sidewall. For further details about 



the numerical code we refer to Verzicco & Orlandi ( 1996 1 and Verzicco & Camussi ( 1997 



2003). In all simulations we calculate the azimuthal averages of the three velocity com- 
ponents and the temperature. Furthermore, in the simulations at Ra = 2.73 x 10 8 we 
placed 32 azimuthally equally spaced numerical probes that provide simultaneous point- 
wise measurements of the temperature and the three velocity components u r , u z 
in the radial, azimuthal, and vertical directions r, <f> and z, respectively at the heights 
0.25L, 0.50L, and 0.75L and a distances 0.45L from the cylinder axis. The azimuthal 
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(a) (b) 



Figure 4. Experimental results on the magnitude of the different Fourier modes and S m for 
Ra = 1.16 x 10 9 and Pr — 4.38 based on the data of 8 azimuthally equally spaced probes at 
z = 0.50L. a) The energy in the different Fourier modes of the azimuthal wall temperature profile 
normalized by the energy present in all Fourier modes as function of 1 /Ro. b) The corresponding 
relative LSC strength Sm- The two dashed vertical lines in both panels indicate the transitions 
between regimes 1 and II, and between II and III, respectively. 



temperature and vertical-velocity profiles measured by the probes are analyzed in the 
same way as the experimental data. 

Figure [5] shows the magnitude of the different Fourier modes of the azimuthal vertical- 
velocity profiles for Ra = 2.73 x 10 8 and Pr = 6.26 based on the data of 8, 16, and 32 
equally spaced probes. The result in panel (a) agrees well with the result obtained by 
Kunnen et al. (2008a). Furthermore, a comparison of the relative energy in the different 



Fourier modes based on the data of 8, 16, and 32 probes, see figure [5] reveals that 8 
probes are insufficient to capture all flow characteristics. Since the results based on the 
data of 16 and 32 probes are very similar, we assume that 16 probes should be sufficient 
to capture all relevant features of the azimuthal profiles. 

Fortunately, the data of 8 equally spaced thermistors are already sufficient to reliably 
calculate S m for non-rotating RB convection, see Stevens et al. (201061. To confirm 
this observation for the RRB case (and validate the experimental results with 8 probes 
discussed in section [3| we therefore calculate S m based on the data of 8, 16, and 32 
equally spaced probes. The result is given in figure [5jd), which shows that the curves for 
the relative LSC strength almost collapse for the three cases. Moreover, the numerical 
data are in good agreement with the experimental result shown in figure |4j including 
the large decrease of S m at the transition between regimes I and II. The small value 
of the relative LSC strength (S m <C 0.5) suggests the absence of the LSC and that the 
azimuthal temperature profile at midheight becomes random in regime II and III. This 



is confirmed by a three-dimensional visualization of the flow, see figure 5 of Zhong et al. 



(2009), where it is shown that in regime II vertically aligned vortices are the dominant 
feature of the flow. We note that figure [5] shows a continuous decrease of the relative 
LSC strength, which is in disagreement with the measurements of Zhong & Ahlers ( 2010 1 
that show a small increase in the temperature amplitude of the LSC between 1/Ro = 
and 1/Ro ~ 0.41. This difference might be caused by the slightly different monitoring 
locations (inside the sidewall in the experiments and at the radial position 0.45L in the 
numerical simulations). 
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Figure 5. Results from DNS on the magnitude of the different Fourier modes and S m for 
Ra — 2.73 x 10 s and Pr = 6.26 based on the data of the numerical probes at z — 0.50L. Panels 
(a), (b), and (c) show the energy in the different Fourier modes of the azimuthal temperature 
profile normalized by the energy present in all Fourier modes based on the data of 8, 16, and 
32 equally spaced probes, respectively. The black, red, blue, and dark green lines indicate the 
energy in the first, second, third, and the additional modes, respectively. Note that when the 
data of 8 probes are used a lot of information about the higher modes is lost. Panel d shows 
the corresponding relative LSC strength based on the data of 8 (red), 16 (blue), and 32 (black) 
equally spaced probes. The two dashed vertical lines indicate the transitions between regimes I 
and II, and between II and III, respectively. 



5. The role of Ekman and Stewartson boundary layers 

In this section we will first discuss the flow characteristics for regimes I and II and 
sketch the influence of background rotation on the mean flow in general terms, and 
the role of the Ekman and Stewartson boundary layers is emphasized. Subsequently, we 
provide physical explanations for the structure of the mean flow as observed for regime 
II. Finally, we discuss the implications for the mean temperature gradient at the sidewall 
of the convection cell. 

5.1. Description of the flow characteristics 

In figure [6] we show the azimuthally averaged temperature and velocity components u z , 
Ut/> and u r for three typical rotation regimes: (the top row) l/Ro — (no rotation, 
regime I), (the middle row) 1/Ro — 0.35 (weak rotation, regime I) and (the bottom 
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row) 1/Ro = 2.78 (moderate rotation, regime II) at Ra — 1 x 10 9 and Pr = 6.4. The 
left-hand side of each picture (r = 0) represents the position of the cylinder axis (dash- 
dotted line), while the right-hand side (r/L = 0.5) corresponds with the sidewall. In 
addition to the azimuthal averaginging, they have also been averaged in time for more 
than one hundred large-eddy turnover times to find the mean circulation hidden under 
the turbulent fluctuations. A close-up of the thermal and flow structure in the bottom 
corner of the tank for the case 1/Ro = 2.78 (see figure |6j:) is shown in figure [7] We note 
that the averaged profiles for lower Ra values are similar to the ones presented in figure 

m 

For the non-rotating case (1/Ro = 0), shown in [6^ a), we observe the signature of the 
LSC. When looking at a vertical cross-section of the domain, the cross-sectional plane 
aligned with the LSC, the mean circulation has an elliptic shape with its major axis 
oriented at some angle with the cylinder axis (see the sketch in figure [8j and also in |Qiu| 



& Tong (2001 1). It is due to this mean tilt of the elliptic LSC, combined with the rotation 
sense as shown in figure [8| that the azimuthal average of this flow is not zero; especially 
upward and downward motions are separated in these azimuthally averaged plots. In the 
top half upward motion (positive u z ) is found near the cylinder axis, while the downward 
flow (negative u z ) is strongest near the sidewall. The opposite situation is found for u z 
in the bottom half. The azimuthal velocity does not show a well-defined mean profile, 
just some small fluctuations. The averaged radial velocity u r has maximal values near 
the intersections of the horizontal plates with the sidewall, consistent with the averaged 
vertical velocity u z near these regions. Note also the weak radially inward flow near the 
sidewall around midheight. 

When a small rotation is introduced [1/Ro = 0.35, figure [6jb)] , a first observation 
is that the mean velocities become smaller (except the averaged azimuthal velocity u$ 
which grows due to organization of the azimuthal flow). The vertical velocity u z still 
has the dominant structuring due to the LSC, but especially near the cylinder axis some 
disturbances appear. For the azimuthal velocity a well-defined structure emerges due to 
the azimuthal averaging. The Coriolis force, deviating the parts of the LSC with mean 
horizontal velocity to the right with respect to the local direction of the (mean) horizontal 
flow, is driving the organization of an average azimuthal flow. The blue regions near the 
corners in the plot of the azimuthal velocity u$ represent the anticyclonic motion induced 
by the mean horizontal outward flow of the LSC (near the bottom and top plates, see 
figure [8]). The sketch of the LSC motion shows that the up (down) going plumes are 
first traveling straight up (down), before the radial inward flow of the LSC sets in. The 
cyclonic motion that is observed in the central sidewall region of the azimuthally averaged 
azimuthal velocity is due to the radial inward flow of the LSC in the top (bottom) half of 
the cell for the upward (downward) branch of the LSC, see figure [8j This radial inward 
flow of the LSC results in a spin-up, i.e. in cyclonic azimuthal motion, by the action of 
the Coriolis force. This mechanism is equivalent to conservation of angular momentum: 
radial inward motion (u r < 0) of fluid parcels results in an increase of u^, i.e. in cyclonic 
motion. Because this radial inward flow is in the top (bottom) part of the sample for the 
up (down) going plumes this results, after azimuthal averaging, in one central sidewall 
region in which a cyclonic motion is found. This means that this particular snapshot, 
cannot be used for interpretation of the movement of individual plumes traveling from 
bottom to top and the other way around. However, in this very schematic picture it does 
seem to be consistent with experimental measurements of the anticyclonic motion of the 



LSC (Hart et al. 2002 Brown & Ahlers 2006a Kunnen et al. 2008a Zhong & Ahlers 



2010) as both upgoing and downgoing plumes get anticyclonic deflection by the mean 
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anticyclonic azimuthal flow near the plates and any effect of mean cyclonic flow becomes 
effective only for z > 2/3H (for upgoing plumes) and z < H/3 (for downgoing plumes). 

At much higher rotation rates, i.e. for larger 1/Ro- values, the LSC no longer exists, 
and the secondary circulation takes on a different appearance. It is obvious from the 
velocity plots that relatively large mean velocities occur in thin regions near the solid 
container walls. A division of the domain into several regions is appropriate, being the 
bulk interior domain, the Ekman layers at the horizontal plates, and the Stewartson layer 
at the sidewalk The Ekman layers have a (non-dimensional) thickness Se — Ek 1 / 2 (with 
the Ekman number defined as Ek = v/(HL 2 ), representing an inverse Reynolds number 
based on the length scale L and velocity scale Lfi). The Stewartson layer at the sidewall 
has a sandwich structure consisting of a thicker outer layer of non-dimensional thickness 



Ek 1 / 4 and a thinner inner layer of non-dimensional thickness 6 3,1/3 = Ek 1 / 3 (see 
e.g. |Stewartson||1957| |1966| |Moore fc Saffman||1969| |van Heijst||l983[ |1984[ |1986[). The 



°S,l/4 



dashed lines in figure |6jc) and figure |7] indicate these typical layer thicknesses. In figure 
[7] the bottom right-hand corner is magnified for the temperature (top panel) and each of 
the velocity components, so that the various boundary layers are more easily recognized. 
Note that since the secondary circulation associated with the Ekman and Stewartson 
layers is weak, the mean velocities u r and u z in the bulk are very small and therefore 
hardly noticeable, as is confirmed by the simulation results. It should be emphasized that 
the secondary circulation found here is considerably different from that found by |Hart"<%£ 



Olsen ( 1999 ) or that of Homsy & Hudson ( 1969 ), both driven by centrifugal buoyancy. In 



both these works the up-down symmetry is broken due to centrifugal accelerations having 
opposite directions near the plates, viz. cold fluid near the top plate accelerates radially 
outward while hot fluid near the bottom plate goes radially inward. Here we disregard 
centrifugal acceleration and up-down symmetry is preserved. It must also be stated that 
the aforementioned circulation is only of secondary magnitude. The still-turbulent field 
of the vortical plumes is dominant. 

The most remarkable change is observed in the distribution of the azimuthal velocity 
u^,: in the bulk of the domain an anticyclonic flow (u^ < 0, indicated by blue in figure 
|6j;) is present, while a thin band of cyclonic flow (u^ > 0, indicated by red) is visible only 
close to the sidewalk Simulations for the parameter range 0.35 < l/Ro < 2.78 (not shown 
here) revealed that for increasing rotation rates (increasing 1 jRo- values) the cyclonic flow 
region becomes confined to a region near the sidewall of increasingly smaller size. These 
observations suggest that in this i?o-regime the Ekman and Stewartson layers play an 
essential role in the mean secondary circulation, even though this circulation is relatively 
weak and even continuously disturbed by the non-steady thermally driven turbulent 
plumes that are present throughout the flow domain. 

5.2. Analysis secondary circulation regime II 

The feature of the weak anticyclonic swirl flow in the bulk of the domain in combination 
with a region of cyclonic swirling motion near the cylinder wall (the u^-plot in figure 
|6jc)) is remarkable. Although the exact mechanism causing this bulk motion is an open 
question, we offer a possible explanation as follows. The heated fluid rising upward from 
the bottom plate and the cooled fluid flowing downward from the top plate result in a 
radially outward mean flow in the bulk of the domain, as sketched in the (r, z) plane 
in figure [9^a) . Although this radial motion is relatively weak, conservation of angular 
momentum (~ rV , with V — fir + the absolute azimuthal velocity) implies that the 
absolute swirl velocity of fluid parcels will decrease, i.e. in the co-rotating frame this fluid 
will acquire an anticyclonic azimuthal motion (u^ < 0). This radial outflow in the bulk 
of the domain is of such small magnitude that it leaves no strong trace in the azimuthally 
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Figure 6. The figures from left to right indicate temporal and azimuthal averages of temperature 
and the velocity components u z , u$, and u r , respectively, for Ra = 1.00 x 10 9 and Pr = 6.4. 
The rows from top to bottom indicate the results for (a) 1/Ro — 0, (b) 1/Ro — 0.35, and (c) 
1/Ro = 2.78, respectively. The left-hand side of each picture is the cylinder axis r = 0; the 
right-hand side corresponds to the sidewall r/L = 0.5. The dashed lines in the pictures of the 
bottom row indicate typical boundary-layer thicknesses: 8e = Ek 1 ^ 2 near the bottom and top 
plates, and 63,1/3 = Ek 1 ^ 3 and 63,1/4 = Ek 1 ^ 4 near the sidewall (83,1/3 is closest to the sidewall). 
Note that the azimuthally averaged azimuthal velocity for 1/Ro = (a) is much smaller than 
for 1/Ro = 0.35 (b) and 1/Ro = 2.78 (c). 



averaged radial velocity u r . Furthermore, the plots of averaged vertical velocity u z show 
no mean motion directed away from the horizontal plates: this transport is localised in 
the small thermal vortices that vanish after azimuthal averaging. 

This picture of the flow is not complete, however, because a second mechanism is active 
simultaneously in the bulk of the domain: the spin-up process (iGreenspan & Howard 



1963). This second process is driven by the Ekman layers (figure ^p): the bottom Ekman 
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Figure 7. Close-up of the bottom-right corner of the figures shown in figurepTc) (1/Ro = 2.78). 
The area displayed is 0.15 < r/L < 0.5; < z/L < 0.1. Color coding ana dashed lines as in 
figure |6|c). The panels from top to bottom indicate the temporal and azimuthal averages of 
temperature and the velocity components u z , u<j>, and u r , respectively. 

layer imposes a suction velocity on the interior given by u z = ^Eku>i in dimensionless 
form, with ljj the vertical component of the rotation of the relative interior flow. To 
clarify, as w/ is anticyclonic, the corresponding vertical velocity according to the suction 
condition is directed from the bulk into the Ekman layer. A similar suction velocity is 
imposed on the interior flow by the upper Ekman layer, also directed from the bulk to 
the boundary layer. There is thus a mean flow from the bulk into the Ekman layers. By 
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Figure 8. Schematic side view of the tilted LSC in the cylinder in the non-rotating or weakly 
rotating (regime I) case. The dash-dotted line is the axis of the cylinder. 




(a) (b) 

Figure 9. (a) Fluid motion due to thermal forcing. The heated (cooled) fluid rising upward 
(flowing downwards) from the bottom (top) plate result in a radially outward mean flow in the 
bulk of the domain, (b) Fluid motion due to the suction of the Ekman layers (spin up). The 
Ekman layers impose suction towards the plates. Mass conservation implies a radial outflow 
inside these layers. The dash-dotted line in both panels indicates the center of the cell (which 
coincides with the rotation axis). 

conservation of mass, a radial outflow in the Ekman layers is needed to carry away the 
flux from the bulk into the Ekman layers. The radial outflow is easily identified in the 
averaged radial velocity plots due to the tiny radial cross-sectional area of the Ekman 
layer, leading to a relatively large mean velocity. 

The radially outward transport through the Ekman layers is returned via a Stewartson 
boundary layer at the cylinder wall. This layer has a small mean radial inward transport 
into the interior distributed along its vertical extent. An indication of this radially inward 
transport is found in the u r plot of figure[7]at (r, z) ss (0.4, 0.025). There is also an internal 
recirculation in this layer: Again, according to the principle of conservation of angular 
momentum rV, with V = flr + u^, the radial inward motion from the sidewall boundary 
layer in the bulk implies a change in the swirl velocity, which is reduced in magnitude 
and approaches zero. 

It can be shown that the matching of the azimuthal velocity component u^(r = 0.5L) to 
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Figure 10. The azimuthally averaged velocity profiles at z = 0.25L for Ra = 1.00 x 10 9 and 
Pr = 6.4 for 1/Ro = 2.78. The azimuthal, radial, and axial velocity components are indicated 
in red, black, and blue, respectively. 



the sidewall requires a Stewartson layer of typical thickness Ek 1 / 4 , in which this velocity 
is given by = — j^u,p(r — 0.5L) exp(^y / 2/L), with L the height of the cylinder, and £ = 
(r — Q.hL^Ek" 1 ^ 4 the stretched radial boundary-layer coordinate (Greenspan & Howard 



|1963 ). When referring to this principal velocity component as being 0(1), the vertical and 
radial velocity components in this layer are much smaller, viz. OiEk 1 ! 4 ) and OiEk 1 / 2 ), 
respectively (in non-dimensional terms). The radial velocity is also correctly matched to 
the sidewall. However, matching of the vertical velocity to the cylinder wall requires the 
presence of a thinner Ek 1 ^ 3 layer inside the Ek 1 / 4 layer. The vertical matching implies a 
principal solution with a vertical velocity u z ~ 0(Ek 1 ^ 4 ) within this layer. However, the 
entire velocity field in the Ek 1 ^ 3 layer carries no net vertical flux; it is just an internal 
recirculation ( Greenspan k. Howard||1963 |. 

We will now describe the internal recirculation in more detail. Where the Stewartson 



Ek 1 / 3 layer meets the Ekman layers, the vertical velocity has a singular structure (Moore 
|fc Saffman|)l969[ |van Heijst|[l983[ |1986[ ): u z ~ ±S(rj) at z = and L, with rj = (r - 
0.5L)Ek~ l/a the stretched radial coordinate in this Stewartson layer, and S the Dirac 
delta function. Signs of these singular eruptions of the Ekman layer fluxes at z = and 
z = L into the thinner Stewartson layer on the sidewall are clearly visible in the u z -plot 
in figure [6|c) , and for the case of the bottom corner in figure [7] Note that these vertical 
fluxes are positive and negative near the bottom and top Ekman layers, respectively. 
Weaker vertical velocities of opposite signs are observed in the thicker Ek 1 ^ 4 layer, see 
figure [(^c). This oppositely signed vertical velocity next to the peak corresponds with the 
internal recirculation of the singular upward flux. Furthermore, the radial branches of 
the internal recirculation can be recognised in the ug plots. Close to the bottom and top 
plates there is strong radially outward flow in the thin Ekman boundary layers, which 
corresponds to strongly negative azimuthal velocity (corner regions near the sidewall of 
figures ^c) and [7]) . Along the sidewall away from the plates there is radial inflow to 
close this recirculation, which, due to the much larger vertical extent, is of such small 
magnitude that it is not found in the plot of mean radial velocity. However, the radial 
inflow leaves its signature on the mean azimuthal velocity: radial inflow corresponds to 
positive ug, which is observed as the orange band near the sidewall in the ug plot of 
figure ^c) . 

To further illustrate the secondary circulation we consider radial mean profiles of ve- 
locity. Figure 10 a) shows the azimuthally averaged velocity profiles u r (r), u^(r), and 
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Figure 11. Boundary layer thicknesses for Ra = 1.00 x 10 9 and Pr = 6.4. The dimensionless 
viscous boundary layer thickness near the plates is indicated by the circles. The dimensionless 
boundary layer thickness of the viscous boundary layers near the sidewall, i.e. the Stewartson 
layers, is indicated by the squares a nd follows the scaling o f the inner Stewartson boundary layer 
thickness, i.e. 63,1/3- Adapted from Kunnen et al.\ |201Qq[ ). 



u z (r) at z = 0.251/ according to a numerical simulation for 1/Ro = 2.78. The radial 
u^-distribution clearly shows the presence of anticyclonic swirl in the interior domain, 
while significant positive swirl exists in the region close to the cylinder wall. The vertical 
velocity shows the peaked structure close to the cylinder wall, indicating the singular 
eruption of the Ekman flux into the Ek 1 ^ 3 layer, with accompanying negative velocity 
directly adjacent to close the recirculation. Furthermore, the radial inflow in the side- 
wall region can be observed here, corresponding to positive ug. In this normalised plot 
this radial-velocity structure may be recognised, yet in absolute sense the radial velocity 
component is much smaller than the vertical and azimuthal contributions. 

The description of the secondary circulation is only qualitative, as a full quantita- 
tive treatment would lead too far for the purpose of the current paper. A quantitative 
treatment of the secondary circulation will be presented in a forthcoming paper. 

Additional evidence that the mean velocity profiles close to the sidewall are indeed 
governed by Stewartson boundary layer dynamics is provided in figure [TT] where it is 
shown that the thickness of the sidewall boundary layers observed in the simulations 



follows the theoretical predictions (Kunnen et al. 2010a 



The role of Ekman and Stewartson boundary layer dynamics on the average mean 
flow is here illustrated with numerical simulations with Ra = 1 x 10 9 and Pr = 6.4. 
Separate sets of simulations with Ra — 2.73 x 10 8 with Pr — 6.26 and Ra = 1 x 10 s with 
Pr = 6.4 essentially revealed the same picture. The explanation provided in this section is 
therefore applicable to a range of Rayleigh numbers. Further studies are needed to explore 
the regime with Ra > lx 10 9 and the role of the Prandtl number. We have also explored 
the role of the Rossby number, in particular by increasing the value of 1/Ro such that 
the flow is dominated by strong rotation (close to and even in regime III). The vertical 
vortex plumes then become stronger and data from azimuthally averaged components 
of velocity reveal an increasing number of alternating, vertically almost homogeneous 
regions of < and > 0. 
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5.3. Influence of the secondary circulation on sidewall temperature measurements 

The above observations of the flow structure enable us to explain the sidewall tempera- 
ture measurements. In figure [12] we have plotted radial distributions of the azimuthally 
averaged temperature and vertical velocity at z = 0.25L. Three main observations from 
the data extracted at this particular level are: a strong vertical mean flow near the side- 
wall (and for the highest rotation rate clearly within the Ek 1 / 3 layer), an increase of 
the mean temperature in the bulk with the rotation rate, and an increased mean tem- 
perature near the sidewall compared to the bulk, regardless of the existence of the LSC. 
We start with the latter observation. In regime I (weak rotation) where the LSC is the 
dominant feature of the flow, the temperature gradient at the sidewall is caused by the 
LSC, which carries warm fluid upwards along the sidewall and cold fluid down in the 
middle (when considered in an azimuthally-averaged plot); and vice versa in the top half 
of the cylindrical domain. In regime II (moderate rotation) the LSC has disappeared 
and vertical vortices form the dominant feature of the flow. In this regime the secondary 
circulation described above causes a flow directed vertically away fro m th e plate close to 



the sidewall, i.e. the singular eruption into the Ek 1 / 3 layer, see figure 12 Id). This means 
that in regime II the secondary circulation carries warm fluid upwards along the sidewall 
in the lower half of the cylinder and cold fluid downwards along the sidewall in the top 
half. The vertical temperature gradient on the sidewall is increased. The enhanced mean 
temperature in the bulk with increasing rotation rate indicates that the vertical tem- 
perature gradient in the bulk increases also with increasing 1/Ro, see for example figure 



15 of Kunnen et al. (2010a). The bulk temperature gradient is caused by the merger of 
vertical plumes (see, e.g., Julien et al. | 1996 b 1999 Legg et aL|2001 Sprague et aL|2006" 



Boubnov fc Golitsyn|1986||Zhong et aL|1993flEcke fc Liu||1998||Sprague et a/.||2006[ )7Te 



the enhanced horizontal mixing of the temperature anomaly of the plumes results in a 
mean temperature gradient. In regime III this temperature gradient in the bulk becomes 
stronger than in regime II and eventually becomes equal to the temperature gradient at 
the sidewall. 



6. Conclusions 

Based on the experimental data and results obtained from DNS we studied the char- 
acteristics of the azimuthal temperature profiles at the sidewall. We find that in regime I 
(weak rotation) the LSC is the dominant feature of the flow. In the sidewall temperature 
measurements this is identified by the strong presence of the first Fourier mode in the 
signal. When the rotation is increased a sudden onset in the heat transport enhance- 
ment is found that indicates the beginning of regime II, the moderate rotation regime. 
In this regime the LSC is replaced by vertically aligned vortices and due to their random 
position the magnitude of the first Fourier mode is very weak, which confirms that the 
LSC has disappeared. This feature is also observed in regime III, where the heat transfer 
decreases, because the vertical velocity fluctuations are suppressed by the rotation. 

When the LSC is the dominant feature of the flow the vertical temperature gradient at 
the sidewall is mainly due to plumes that travel close to the sidewall because they travel 
with the LSC. However, also in regime II and III, in which the LSC is absent, there 
is a strong temperature gradient at the sidewall. In these regimes, where the vertical 
vortices are the dominant feature of the flow, the observed vertical temperature gradient 
along the sidewall is mainly due to the secondary flow. The secondary flow, driven by 
the Ekman boundary layers near the plates, causes a recirculation in the Stewartson 
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Figure 12. (a) Azimutally averaged temperature and (b) vertical velocity profiles as function of 
the radial position at the height z = 0.25L for Ra — 1.00 x 10 9 and Pr = 6.4. Black, blue, and 
red indicate the data for 1/Ro — 0, 1/Ro = 0.35, and 1/Ro = 2.78, respectively. The dashed 
lines indicate typical Stewartson boundary layer thicknesses (<5s,i/3 = Ek 1 ^ 3 and 83,1/4, = Ek 1 ^ 4 
near the sidewall for 1/Ro = 2.78 (63,1/3 is closest to the sidewall). Note that for 1/Ro — 2.78 
there is a strong upward flow of fluid in the inner Stewartson Ek 1 ^ 3 layer, which causes a 
vertical temperature gradient at the sidewall. In panel (a) the increasing mean temperature in 
the bulk with increasing rotation rate indicates that the vertical temperature gradient in the 
bulk increases with increasing 1 / Ro. The temperature gradient in the bulk for this case is shown 
in figure 15 of Kunnen et al. (2010a I. 



boundary layer on the sidewall with upward (downward) transport of hot (cold) fluid 
close to the sidewall in the bottom (top) part of the cell. 

It is remarkable that the secondary flow in turbulent RB convection, which is observed 
after time-averaging of the flow field, is so well described by linear Ekman and Stewartson 
boundary layer theory. It seems that the mean laminar flow profiles prevail, hidden under 
highly turbulent fluctuations. Thus the knowledge about laminar flow dynamics can still 
be highly relevant in the study of turbulent flows. An example is the Grossmann-Lohse 
(GL) theory 



see 



Ahlcrs et al. 2009 



for an overview), in which the kinetic energy and 
thermal variance dissipation rates have been decomposed into boundary-layer and bulk 
contributions. Scaling- wise and in a time-averaged sense a laminar Prandtl-Blasius profile 
is assumed close to the horizontal plates. The GL theory successfully describes the Nu 



and Re number dependences on Ra and Pr. Recently, Zhou & Xia (2010) and Zhou 



et al. (20101 have shown that the laminar PB profile is indeed hidden under the turbulent 



fluctuations. 
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